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Abstract 

m 



We suggest that supersymmetric dark matter be used to set the energy scale of a linear e + e 



collider. Assuming that the lightest supersymmetric particle (LSP) is a stable neutralino \i 
as in many incarnations of the MSSM with conserved R parity, previous calculations that 
include coannihilation effects have delineated the region of the (wi/2,?no) plane where the 
LSP cosmological relic density lies in the preferred range 0.1 < Q x h 2 < 0.3. We evaluate 
here the total cross section for e + e~ — > visible pairs of supersymmetric particles, for different 
values of mi/2 and mo, and investigate how much of the dark matter region can be explored 
by e + e~ colliders with different centre-of-mass energies Eqm- We find that a collider with 
Ecm — 500 GeV or 1 TeV can only explore part of the cosmological region, and that a collider 
with Ecm = 1.5 TeV with sufficient luminosity can explore all of the supersymmetric dark 
matter region. 
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1 Introduction 



One of the most promising ideas for a high-energy accelerator to complement the LHC is 
a linear e + e~ collider (LC) with a centre-of-mass energy Eqm in the TeV range fl[]. The 
crucial parameters of such a LC are Ecm and the luminosity. The optimal choice of Ecu 
is constrained by technology and cost, but should be driven by physics arguments based 
on the accessibility of physics thresholds. One established threshold in the energy range of 
interest is that for e + e~ — > it at about 350 GeV @. A second threshold likely to be in this 
energy range is that for Higgs boson (H) production via the reaction e + e~ — > H + Z ||. 
For some years ||], the precision electroweak data have favoured a relatively light Higgs 
boson, as suggested independently by supersymmetry. The most recent indication is that 
Mh 200 GeV at the 95 % confidence level || , corresponding to a H + Z threshold below 
about 300 GeV. 

Since supersymmetry is widely considered to be one of the most promising possible low- 
energy extensions of the Standard Model, it is desirable that any new collider offer good 
prospects of detecting at least some supersymmetric particles, as is the case of the LHC 0. 
The physics argument that has usually been employed to estimate the sparticle mass scale rh 
has been that of the naturalness of the gauge hierarchy, which suggests that m < 1 TeV [[7]. 
A supporting argument has been the concordance of the gauge couplings measured at LEP 
and elsewhere with the predictions of supersymmetric Grand Unified Theories (GUTs) |§. 
However, this argument is sensitive only logarithmically to rh, and is also vulnerable to 
GUT threshold effects due to particles beyond the minimal supersymmetric extension of the 
Standard Model (MSSM). The agreement of the Higgs mass range favoured by the precision 
electroweak data with that calculated in the MSSM U is also encouraging, but is again only 
logarithmically sensitive to rh, and hence unable to specify it with any accuracy. 

An independent argument for new physics around the TeV scale is provided by calcula- 
tions of cold dark matter, which yield naturally a freeze-out density in the cosmologically 
allowed range: VtcDNih 2 < 0.3 (where Q = p/p c , the critical density, and h is the Hubble 
expansion rate in units of 100 km/s/Mpc), and that preferred by theories of structure for- 
mation: 0.1 < QcDMh 2 , if the mass of the cold dark matter particle is < 10 TeV |[10|| . The 
upper limit on flcDAih 2 is fixed by the age of the Universe. For Q tot < 1, a lower limit on 
the age of the Universe of 12 Gyr implies an upper limit fl m h 2 < 0.3 on the total matter 
density, and hence VLcdm < Q m - This argument does not rely on the high-redshift supernova 
observations ||TTJ, but they do support it. 

A serendipitous prediction of flcDM is provided by the MSSM with R parity conserva- 
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tion [12J, if the lightest supersymmetric particle (LSP) is the lightest neutralino Xi as i n 



many versions of the MSSM. Indeed, it has been shown [J3| that the most 'natural' choices 
of MSSM parameters, from the point of view of the gauge hierarchy, yield a relic LSP den- 
sity in the astrophysical and cosmological region 0.1 < Q x h 2 < 0.3. In this case, detailed 
calculations of the relic LSP abundance yield flcDM < 0.3 only for m x < 600 GeV |TJ]]. An 



essential role in this relic density calculation is played by x ~ & coannihilation effects when 
the LSP is mainly a gaugino, which increase significantly the upper limit on the LSP mass 



quoted previously |]15 | 



The idea we propose in this paper is that the relic density calculation be used to specify 
the likelihood that a LC with given Ecu will be above the sparticle pair-production thresh- 
old, and able to detect at least some supersymmetric cross section. The answer is necessarily 
higher than Ecu = 2 x m™ ax , since the process e + e~ — > XX is n °t directly observable in 
models with a stable neutralino LSP x- On the other hand, as we discuss in more detail 
below, m™ a:r ~ 600 GeV is attained when m x ~ rrif, with m^,mg not much heavier, so one 
might expect that a LC with Eqm not far above 1200 GeV should be sufficient. As we show 
in more detail below, a LC with Ecu = 500 GeV or 1 TeV would only be able to detect 
supersymmetry in a fraction of the preferred dark matter region of MSSM parameter space. 
A LC with Eqm — 1-5 TeV would probably cover the preferred region, but might miss some 
part of the x ~ & coannihilation 'tail' at large mi/2, depending on the luminosity it attains. 
A LC with Eqm = 2 TeV would, on the other hand, be able to cover all the cosmological 
region with a comfortable safety margin in terms of cross section, kinematic acceptance and 
astrophysical uncertainties. 

2 Summary of LSP Density Calculations 

We assume R parity is conserved, otherwise there would be no stable supersymmetric 
dark matter to interest us. We work within the constrained MSSM, in which all the 
supersymmetry-breaking soft scalar masses are assumed to be universal at the GUT scale 
with a common value mo, and the gaugino masses are likewise assumed to be universal with 



common value mi/2 at the GUT scale |16j. The constrained MSSM parameters are chosen so 
as to yield a consistent electroweak vacuum with a value of tan j3 that is left free. The LEP 
lower limits on MSSM particles, including the lightest Higgs boson, suggest that tan/3 > 3, 
so we consider this and the higher value tan/5 = 10. We consider two possible values of the 
trilinear soft supersymmetry-breaking parameter: A = 0, -m^, the latter being the value 
for which the constraint that the lowest-energy state not break charge and colour (CCB) is 
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weakest [0 , consistent with parameter choices out to the point at the tip of the cosmological 
region. 

When calculating the relic density of LSPs Xi h is assumed that they were in thermal 
equilibrium prior to freeze-out at some temperature Tf. The relic density after freeze-out 
is then determined by the competition between the expansion rate of the Universe and the 
neutralino annihilation rate. Ultimately, the relic density is inversely related to the effective 
annihilation cross section <J e ff, which falls off as the square of the supersymmetry breaking 
scale. Thus, as the supersymmetry breaking scale is increased, the annihilation cross section 
decreases and the relic density increases. This is why an upper limit to the relic density 
puts an upper limit on the sparticle mass scale, and on the mass of the neutralino LSP, in 
particular. In regions where the neutralino is mainly a gaugino (usually a bino), as in many 
models of interest, such as those with GUT-scale universality relations among the sparticle 
masses, the annihilation rate is dominated by sfermion exchange. As one approaches the 
upper limit on the neutralino mass, the cross section is maximized by taking sfermion masses 
as small as possible: in this case, the sleptons i are nearly degenerate with the neutralino 
LSP X Q 

When the LSP is nearly degenerate with the next-to-lightest super symmetric particle 
(NLSP), it is known |TB[ that new important coannihilation channels must be included to 
determine the relic neutralino density. Thus, in addition to the self-annihilation process 
XX ~^ anything, the effective annihilation cross section includes important contributions 
from coannihilation processes involving slightly heavier sparticles X, Y: %X — > anything, 
XY — > anything, weighted by the corresponding Boltzmann density suppression factors: 



eff ~ <t(xx) + E ie - (m ^^ )/T %(xX) + Z x ye- (m x +m Y- 2m *V T fcT(XY) (1) 



In the parameter region of interest after taking into account the LEP exclusions of light 
sparticles, the most important coannihilation processes are those involving the NLSP f and 
other sleptons: e, /2, which are all taken into account in the following analysis |14j]. Several 



of these coannihilation cross sections are much larger than that for XX annihilation close to 
threshold, because they do not exhibit P-wave suppressions. Therefore, coannihilation is an 
essential complication. 

As noted above, since the resulting LSP relic density Q x h 2 increases as er e // decreases, and 
since a e ff decreases as m , m 1( / 2 increase, one expects generically that Q x h 2 should increase 
with increasing m , ^1/2- This simple correlation is complicated in the presence of nearby 
s-channel Z° and Higgs poles in the annihilation cross sections, but the LEP exclusions 
lr The GUT universality conditions then imply that the squarks are considerably heavier. 
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now essentially rule out this possibility |I9| . As mentioned earlier in the paper, the preferred 
range of cold dark matter density is 0.1 < flcDMh 2 < 0.3. It is possible that all the cold dark 
matter may not consist of LSPs x> so we can a t best assume that Q x h 2 < QcDAdh 2 < 0.3. 
However, this upper limit on tt x h 2 is sufficient to infer an upper limit on m 0l mi/2 f\. In 
TA\, the values of the two key supersymmetry-breaking inputs m^mxii were constrained 



so that neutralino relic density should fall within the desired range. Roughly speaking, 
when mi/2 & 400 GeV, there is a relatively broad allowed range for m between about 
50 and 150 GeV, depending on tan/3,v4 and the sign of fi. For values of m 1 / 2 > 400 GeV, 
coannihilation becomes important, and mo is restricted to a relatively narrow range of typical 
thickness 5rrio ~ 20 GeV. The maximum value of mi/2 is determined by the point where 
there is no longer any value of m , such that the neutralino mass is less than the tr mass 
and QcDMh 2 < 0.3. This occurs when mi/ 2 — 1400 GeV, corresponding to the neutralino 
mass of about 600 GeV mentioned previously. 

This is the essence of our argument that the relic density calculation can be used to specify 
the e + e~ collider energy required to produce sparticles. 

The upper limit to the neutralino mass including coannihilation effects of m x < 600 GeV 



is relatively insensitive to such MSSM parameters as tan/3 and A. As in [lj|, we consider 
here the two cases tan (3 = 3, 10, and initially set A close to the weak-CCB value A = -m 1( / 2 . 
As mentioned earlier, the upper limit on m x implies that the threshold for pair-producing 
sparticles must be at least Ecu = 1200 GeV. In fact, when the limit m x ~ 600 GeV is 
reached, one also has m x = m^, where the NLSP f\ is the lighter stau mass eigenstate, 
so the threshold for the reaction e + e~ — ► f + f~ is also ~ 1200 GeV. Moreover, the mass 
of the e.R is also not far above 600 GeV, so the threshold for e + e~ — > e^e^ is also not far 
beyond 1200 GeV. In addition, it is easy to check that even if one allows m x < m? 1 , which is 
possible if m x < 600 GeV, the threshold for e + e~ — > ? + ?~ is never above 1200 GeV. These 
arguments are all suggestive that Ecm — 1200 GeV may be sufficient for an e + e~ linear 
collider to observe supersymmetry, but any such conclusion must hinge upon the analysis of 
the observability of the sparticle pair-production cross section that we undertake next. 



2 On the other hand, the lower bound on flcDMh 2 > 0.1 cannot be transferred to a lower bound on fl x , 
and hence there are no corresponding lower bounds on mo, mi/2, except for those imposed by slepton searches 
and/or the requirement that the f not be the LSP. 
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3 Analysis of Sparticle Pair-Production Cross Sections 



In order to determine the region of the (m ,m^ 2 ) plane that can be explored with a linear 
e + e~ collider of given Ecm, we have calculated the total observable production cross section 
for the pair production of sparticles e + e" — > XY, where X and Y are not necessarily 
particle and antiparticle [p0|| . In this context, 'observable' means that we do not include 
pair production of the LSP: e + e~ — ► XX- Nor do we include sneutrino pair production: 
e + e~ — > vi>, although some v decays might be visible. Also, the production cross sections 
for heavier neutralinos x' i e -g-? e + e~ — > xx'i are corrected for invisible x' decay branching 
ratios. Finally, we assume that the ordinary particles emitted in a sparticle decay chain are 
observable only if the mass difference AM > 3 GeV. 

We assume an integrated luminosity £ = 100 fb -1 In order to estimate the corre- 
sponding sensitivity to the new-physics cross section a, the relevant quantity is B = y/cr bg / e, 
where a^g is the residual cross section for background processes, and e is the signal-detection 
efficiency As usual, a five-standard-deviation discovery is likely if a > 5 x B / \^C, whereas, 
in the absence of any observation, new-physics processes with a > 2 x B/\fC will be ex- 
cluded at about the 95 % confidence level. At LEP 2, for mass differences between the 
produced sparticle and the LSP that are not too small, the background to searches for 
charginos x ± an d sleptons is mainly due to production, and typical values for B were 
in the range 3 4- 6 (fb + ^). At the LC we expect cleaner background conditions for both 
slepton and chargino searches because the W ± should be more easily distinguishable, and 
also cr(e + e~ — > W + W~) is smaller. It is therefore likely that B is smaller than at LEP 2. 
We adopt a conservative approach and scale B roughly by cr(e + e~ — ► W + W~), taking 
B = 2 (fb + 3), which gives a lower limit on the discoverable cross section of 1 fb, and an 
exclusion upper limit of 0.4 fb. 

Fig. 1 shows the physics discovery reach in the (mo,mi/2) plane for tan/5 = 3,10 pro- 
vided by the processes e + e~ — > £ + £~ , neutralinos and charginos x + X~ f° r collisions at e + e~ 
collisions at E C m = 500, 1000, 1250, 1500 GeV, compared with the allowed cold dark matter 
region (shaded). The solid lines in Fig. 1 correspond to the estimated discovery cross section 
of 1 fb for e + e _ — > £ + £~ , and the broken lines to the kinematic limit m x ± = Ecm/2- We see 
no big differences between the plots for the different signs of fx, nor indeed for the different 
values of tan (3. We note that e + e _ — ► £ + £r (solid lines) provides the greatest reach for each 
of the values E C m = 500, 1000, 1250, 1500 GeV studied, and that chargino pair production 
e + e~ — > x + X~ (broken lines) becomes progressively less important as Ecu increases. 

We see in Fig. 1 the extent to which the region favoured by the cosmological requirement 
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Cross section limit a,. = 1 fb 




Figure 1: Discovery sensitivity in the (mo, m^) plane for tan (3 = 3 (left panels) and 10 
(right panels) provided by searches for e + e~ — > £ + £~ and neutralinos (solid lines) and x + X~ 
(broken lines) for collisions at Ecm = 500,1000,1250,1500 GeV. The allowed cold dark 
matter regions are shaded. The top (bottom) panels are for fi < (>)0 ; and the value A = 
—mi/2 is used. 
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that 0.1 < Q x h 2 < 0.3 may be covered by LC searches at different energies. In particular, 
about a half of this region is covered by sparticle searches at Ecm = 500 GeV, a somewhat 
larger fraction (but not all) is covered at Ecm = 1000 GeV, and full coverage of the favoured 
region is approached only when Ecm — 1500 GeV The reason why more than 1200 GeV 
is required is the P-wave threshold suppression for the observable processes with the lowest 
thresholds near the point of the cosmological region, namely the reactions e + e~ — > t^l^. 

Fig. 2 shows as three-dimensional 'mountains' the full observable sparticle cross section 
for tan/3 = 10 and /i > for E C m = 500, 1000, 1250 and 1500 GeV, including also other pair- 
production processes. The irregularities in the outline of the three-dimensional 'mountain' 
plot correspond to the opening up of different sparticle pair-production thresholds. We see 
again that E C m = 500 GeV is not adequate to cover much of the cosmological region, that 
Ecm = 1000 GeV does not cover a significant fraction of the high-mx/2 tail opened up by 
coannihilation, and that Ecm > 1500 GeV covers the cosmological region. We find similar 
features for tan/3 = 10 and fi < 0, and also for tan/3 = 3 and both signs of // (not shown). 

We now return to the tip of the cosmological tail, which occurs when m x ~ 630(610) GeV 
for tan/3 = 3(10) for our default option Q x h 2 < 0.3, and explore in more detail how much 
Ecm beyond 1200 GeV is required to be sure of detecting supersymmetry. Fig. 3 shows the 
contributions to the effective observable cross section from the dominant reactions e + e~ — > 
it,R^I,R- Close to threshold, only pair production of the £r states is accessible, which exhibits 
a P-wave suppression. The associated-production process e + e~ — > tifiR kicks in at somewhat 
higher energies, and rapidly dominates, because of its S'-wave threshold. This is the origin 
of the kink seen in the rise of the total cross section in each of the panels of Fig. 3, where 
the discovery and exclusion sensitivities are also shown as horizontal broken lines. We see 
that Ecm only just above 2m x ~ 1200 GeV is not sufficient for sparticle discovery, because 
of the small observable cross section. We recall that, for our assumed integrated luminosity 
of 100 fb _1 and detector performances, the discovery cross-section limit would be 1 fb, as 
indicated by the upper horizontal broken line in Fig. 3. Of course, this may be altered by 
different assumptions on the integrated luminosity and/or detection efficiency Q. 

Each of the panels in Fig. 3 exhibits alternative curves to be compared with our default 

choices Q x h 2 = 0.3 and A = 0. The curves for Q x h 2 = 0.4 are for instruction only. In this 

case, one finds m x < 740(710) GeV for tan/3 = 3(10), but it is very difficult to reconcile such 

a large value of Q x h 2 with the emerging measurements of cosmological parameters (]. In fact, 

3 We note in passing that a LC with Ecm = 500 GeV would have seemed perfectly adequate if coannihi- 
lation were not taken into account. 

4 We note, in particular, that higher luminosities may be achievable at higher Ecm- 

5 For the record, for Q, x h 2 < 0.5, the upper limit on the neutralino mass increases to m x < 830(800) GeV 
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tanp=10, ji>0, A fi = -m 




Figure 2: The observable sparticle pair-production cross section in the plane for 

tan/? = 10, fi > and A = -m 1/2 , for E CM = 500, 1000, 1250 and 1500 GeV. Note that the 
vertical scale is logarithmic, and that cross-section contours are indicated by changes in the 
shading of the cross-section 'mountain'. The cosmologically-pref erred domain of (mi/2, m o) 
is visible in projection in the top two panels: in the bottom two panels, it is obscured by the 
'mountain'. 
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Figure 3: The total observable cross section for e + e~ — > £ + £~ processes, as a function of 
Ecm, f or the points at the tips of the region of mi/ 2 ,m allowed for Q x h 2 < 0.2,0.3 (our 
preferred choice) and 0.4, with the usual choices tan/3 = 3, 10 and both signs of /i. The solid 
lines are for A = -m^, and the dashed lines for A = 0. The horizontal broken lines are 
our estimates of the possible discovery and exclusion cross sections. 
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we actually believe that allowing Q x h 2 < 0.3 is already quite conservative. For the preferred 
observational value h ~ l/v2, this would correspond to Q x < 0.6, which extends far beyond 
the currently favoured range Q x < 0.4. If instead one enforces Q x h 2 < 0.2, one finds that 
the maximum value of the LSP mass becomes m x ~ 520(500) GeV, for tan/3 = 3(10) and 
Ecm = 1500 TeV would be adequate, as seen in Fig. 3. Indeed, in this case, E C m = 1200 GeV 
would be sufficient to cover all the region of the (mo ; m i/2) favoured by cosmology. We also 
show in Fig. 3 comparisons between the cross sections at the extreme points for A = and 
—mi/2- Our conclusions are clearly insensitive to the ambiguity in the choice of A. 

4 Conclusions 

Finally, we show in Fig. 4 the fraction of the cosmo logically-allowed region of the (mi/2, ttlq) 
plane that can be explored by a LC as a function of the accessible limiting cross section au m , 
for different values of E C m- When the detector perfomances are specified, the values of cr^ m 
correspond to different values of the available luminosity, as indicated. We see in Fig. 4 that 
a LC with Ecm = 1.5 TeV would cover all the cosmo logical region if ou m < 5 fb [], and one 
with Eqm = 1-25 TeV if <Jn m < 0.5 fb. On the other hand, a LC with Ecm = 1 TeV could 
never cover all the cosmological region, and a LC with Ecm = 0.5 TeV covers ~ 60 % of it [} 
The conclusions to be drawn from this analysis are somewhat subjective, since they de- 
pend how much you are prepared to bet at what odds. It could well be that new cosmological 
data might inform better your choice. For example, you could become more sanguine about 
the prospects for a lower-energy LC if the upper limit on Q x h 2 could be decreased to 0.2. Our 
point in this paper has been to establish that there is a phenomenological connection between 
the LC energy and supersymmetric dark matter, and we believe that Fig. 4 summarizes the 
best advice we can offer at the beginning of the third millennium. 
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for tan /3 = 3(10). 

6 A LC with Ecm — 2 TeV would always cover all the cosmological region, even for a very pessimistic 
assumption on ou m . 

7 Fig. 4 is plotted using a linear measure for the cosmological region. The prospects for lower-energy 
machines would seem brighter if one used a logarithmic measure of the parameter space, e.g., using this 
measure, a LC with Ecm — 0.5 TeV would cover over 80 % of the cosmological region. 
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Figure 4: The discoverable fractions of the region of the (m 1 / 2 , m ) plane allowed by cosmology 
for 0.1 < Q x h 2 < 0.3 that is accessible to a LC as a function of accessible cross section o- Um , 
for different values of E C m, for each of our scenarios tan/? = 3, 10 and both signs of \x, 
assuming A = -m^. Also indicated is the correspondence between luminosity and au m for 
the detector performances assumed in this paper. 
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